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We develop an analytical technique to evaluate the magnetic helicity in the debris from sphaleron 
decay. We show that baryon number production leads to left-handed magnetic fields, and that the 
magnetic helicity is conserved at late times. Our analysis explicitly demonstrates the connection 
between sphaleron-mediated cosmic baryogenesis and cosmic magnetogenesis. 



I. INTRODUCTION 



A cosmic phase transition in which magnetic monopoles are produced necessarily produces cosmic magnetic fields, 
as these are sourced by the monopoles. The electroweak model also contains magnetic monopoles [Ij that can source 
magnctic fields However, a quantitative estimate of the magnetic field is difficult to obtain because electroweak 

monopoles are confined to electroweak antimonopoles by Z-strings [1] , and estimates of the monopole density based on 
topological considerations cannot be applied. Fortunately there is a way to sidestep this difficulty because magnetic 
monopole-antimonopole pairs, in the guise of "sphalerons" [6], play a crucial role in the violation of baryon number 
in the electroweak model [5 . Hence monopole-antimonopole pairs are a common factor in baryogenesis and magne- 
togenesis. This fact allows one to estimate the hehcity of the magnetic field by relating it to the baryonic density 
[II [8]. The connection between baryons and magnetic fields may provide an opportunity for tests of fundamental 
particle physics and cosmology at the electroweak epoch via observations of the topological and spectral properties of 
a primordial magnetic field. 

To test the production of magnetic fields during baryon number violating processes, sphaleron decay was studied 
in pl (also see [10]) by numerically evolving the electroweak equations of motion on a lattice. Debris from sphaleron 
decay was analysed for magnetic heficity and the connection between baryon number violation and generation of 
magnetic helicity was confirmed. The aini of the present paper is to analytically study the sphaleron decay process 
and to derive the generation of magnetic hehcity in a transparent manner. 

Ahhough the fuU dynamics of sphaleron decay involves the complete electroweak equations of motion, we are able 
to extract the dynamics relevant to the electromagnetic sector foUowing a scheme similar to Ref. [TT] . Our strategy is 
to first consider the "S'Í7(2) sphaleron" which is an expUcitly known saddle-point sohition in the electroweak model 
with the hypercharge coupHng, g' , set to zero [^. Not only is the SU{2) sphaleron known expUcitly but a path in 
field space that connects the sphaleron solution to the vacuum solution has also been constructed [6]. Let /i be the 
parameter along this path. We wiU consider the decay of the sphaleron along this particular path in field space. 
Hence the time dependence of the fields is given by the (unspecified) function /x(í). With this model of the decaying 
sphaleron, we now turn on the hypercharge coupUng, g' , assuming it to be smaU. Then the decaying sphaleron 
produces electromagnetic currents, which produce electromagnetic fields. We calculate the magnetic heUcity in these 
fields and show analyticaUy that it is conserved at late times and asymptotes to some non-zero value which depends 
on /Li(í). We can also show that the sign of the magnetic heUcity depends on which side of the saddle the sphaleron 
roUs down. TUus the sign of the magnetic heUcity is directly related to whether baryons or antibaryons are produced 
during the decay of the sphaleron. Since we know excess baryons are produced in the universe, cosmic magnetic fields 
carry left-handed magnetic heUcity. Furthermore, the magnitude of the heUcity is directly proportional to the cosmic 
baryon number density as discussed in [3-(Ü]. 



The SU{2) sphaleron is a saddle-point solution to the static electroweak equations with the hypercharge coupUng, 
g' , set to zero. For our purpose, it is necessary to consider a path of field configurations that starts with the sphaleron 
and ends at the vacuum. Such a path was constructed by Manton [B] in order to implement a Morse theory argument 
indicating the existence of a sphaleron solution. The path in field space is now written in spherical coordinates, 
(r, 9, (p), and /i is a parameter along the path 
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where g is the SU{2) coupUng, v/\/2 is the vacuum expectation vahie of the Higgs field, h{r) is the Higgs profile 
function and f(r) the gauge profile function. 
The asymptotic fields are given by 

^^{ii.e,c^)^( sin^sinee'^ \ 
^^' '^^ \^e"'''(cos^ + ism/x cosé'^y 
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The sphaleron solution is obtained if we set ¡i = tt/2 and the vacuum is obtained for fi = or fi = ir. The path from 
/i = 7r/2 to /i = describes the sphaleron rohing down one side of the saddle, and the path from /i = 7r/2 to /í = tt 
describes roU down on the other side of the saddle. The former path describes the creation of antibaryons while the 
latter describes the creation of baryons. 

The dimensionless sphaleron profile functions are found by numericahy solving the static electroweak equations of 
motion. However, the essential features of the profiles are captured by the simple ansatz called "Ansatz o" in [TT], 

II, II r > rf 

'yt^' ■ (6) 
1^1, it r>rh 

The length scales r/ and r^ are related to the vector and scalar masses in the model. For simphcity we wiU take 



rf = r/j, when we numericaUy evalute certain integrals in Sec. VI 



III. ELECTROMAGNETIC CURRENTS 



To modcl the decay of the sphaleron, we take the parameter /i to be a function of time, /i(í), with /i(0) = 7r/2 
and /i(oo) = tt. The Higgs and SU{2) gauge fields are stiU given by Eqs. ([T]). As discussed in [TT], we can now 
consider a small hypercharge coupUng constant, g' , and evaluate tUe electromagnetic field perturbatively in g' . TUen 
the decaying sphaleron wiU be a source for an electromagnetic field. If we denote tUe electromagnetic gauge field by 
A^, MaxweU equations in tUe Lorenz gauge give 

{d'¡ - V^)^^ = , (7) 

wUere J'^ is tUe electromagnetic current produced by tUe decaying spUaleron. To lcading order in tUe UypercUarge 
coupUng constant, g' , tUe electromagnetic current is equal to tUe UypercUarge current given by 

Jl" = -i^ [$t£)M$ _ (£)*' $)t$] . (8) 

To calculate tUe magnetic UeUcity produced by spUaleron decay we only need to find tUe spatial components of 
A^. Hence we only need tUe spatial components of tUe current J'^. Eq. ([s]) gives tUe contravariant current vector 
components; we prefer to work witU tUe "ordinary" components (see Sec. 4.8 of ¡12j ) tUat are obtained by dividing 
by tUe square root of tUe metric factors, namely, r and rsiní? in tUe 9 and </> components. TUe ordinary components 
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of the current are 
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e Je(t,r) (10) 

^^^[/(i-M'cos^M+^'a-/)] 

OMt.r), (11) 

where the prime on h denotes diíferentiation with respect to r. 

The spherical components of the current, Ji, as given in Eq. ([8| only involve spatial derivatives of $ and the 
spacehke components of the gauge fields. Hence the three-current does not involve time derivatives of /i(í) or the 
time-component of the gauge fields. Also, the zeroth (time) component of the current is unspecified since we wiU not 
need it in our calculation. We assume it is such as to ensure current conservation. 

The angular form of the currents is of interest since there is an azimuthal component, J^, and the Jr, Je components 
circulate in the constant (j) plane. Hence the J^ current resembles circular current-carrying wires in planes parallel to 
the a;y— plane and the J^, Je components resemble current-carrying wires wrapped on a toroidal solenoid. The two 
components of the current are hnked and hence the current is hehcal as measured by J • V x J. 

Aheady we can see the relation of the sign of the helicity and the decay path taken by the sphaleron. The path 
from /i = 7r/2 to /i = TT is related to the path from /z = 7r/2to/x = 0by the transformation /i — tt — /t. Under this 
transformation J^ is unchanged but Jr and Je change sign. Thus the linking number of the current components is 
reversed. This impHes that the magnetic helicity produced by the currents wiU also be reversed. 

We wiU eventually need to Fourier transform the current. For this we find it most convenient to work with Cartesian 
components 

Ja;(í,x) = J^siní? cos (f) + J e cos 9 cos — J^ sin (/) (12) 
Jy{t,Tí) = J^siní? sin + Je cos sin^ + J^cos^ (13) 
Jz(í,x) = JrCosO ~ JesvLíO . (14) 

The Fourier transform is implemented by using the plane wave expansion (see Sec. 11.3 of ¡13j ) 

e^k-x ^ V(2n+ l)í" V e„p^^cos(m((/.-t;))P™(cosu)P„"(cos0)j„(fcr) , (15) 
^-^ ^-^ (n + m)! 

n=0 m=0 ^ ' 

where 

X = r^siné* cos (/), siné' sin0, cos 6*) 

k = fc(sinit cos sinw sinu, cosu) , (16) 

and jn{kr) is the spherical Bessel function of order n, P^i-) is the associated Legendre function, and = 1, e„i = 2 
for m > 1. 

Now, for example, the Fourier transform of the x— component of the current is 

J,(í,k) = dn drr^ e''""J,(í,x) . (17) 



In this we insert e*'^'^ from Eq. (15l and the current from (12), using ([9|-(11|. After quite a bit of algebra, we find 

Jx(t,k) = —ATisÍTíu[cos.ucosvK2{t,k)+is.viivKi{t,k)] (18) 

g'v^ 

Jy^t^'k) = — 47rsinu [cosusinw iÍ2(í, ^) — í cosw ifi(í, fc)] (19) 

J,(í,k) = ■l— — [K^{t,k) + {l-?,cos'u)K2{t,k)], (20) 
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where 



Ko = 
Ki = 
K2 = 



and 
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drr'^ji{kr)J^{t,r) 



drr^j2{kr)Jr+e{t,r) 



Jr+e{t,r) = Jr{t,r) + Jg{t,r) 



(21) 
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IV. MAGNETIC HELICITY 



The magnetic hehcity is dehned by 



H{t) = J d^x A-B . 



(23) 



The procedure to evaluate magnetic hchcity is now straightforward though cumbersome. We find A from MaxweU's 
equation, then the magnetic field, B = V x A, and finally thc resuhing magnetic hehcity. In terms of the fields in 
Fourier space (denoted by over tilde's) 



Á(fc) = - 



J(fe) 

W2 



where 3{k) is given in Eqs. (18)-(20l , and 



B(fc) = -ik X A(fc) 



(24) 



(25) 



For us it is more useful to work with the three dimensional (spatial) Fourier transform of the fields while not trans- 
forming their time dependence. This is because thc timc dependence occurs via the unknown function /lí(¿). To go 
from the frequency dependcnt function to the time dependent function, we use 
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27rsin(fc(í- r)) 
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where we have chosen a contour prescription so as to get retarded solutions. 



(26) 



is 



We use Eqs. (24) and (25) in the hehcity integral (|23|) and then integrate over the frequencies using ( |26[ ). The resuh 
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The triple product can be evaluated using the exphcit expressions for the current components in Eqs. (18)-(20) 



Further simphfication involves considerable algebra and integration over trigonometric functions. We also use the 
spherical Bessel function relation 



.h^x) = -ji{x) -jo{x) , 

X 

and perform the integrations over the angular coordinates u, v in fc— space to eventually get 

^W-y (^) [no{t) + ni{t)], 
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/t pt pOO fOO 
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X— / dk[cos{k{T — t')) — cos{k{2t — T — T'))]ji{kr)ji{kr') 



(31) 



The integrations over k are performed using the exphcit expressions for the spherical Bessel functions in terms of 
trigonometric functions 



jo{x) 



sin(a;) 



(32) 



We find 
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and 



where 



with 



sin(.T) cos(x) 
Ji{x) = ^ 



dk k cos{kS)ji{kr)jQ{kr') ^ 
Sr^r' 



Po{S,r,r') 



Po{S, r, r') = r'[So{S, r, r') - So{S, -r, r') - So{S, r, -r') + So{S, -r, -r')] 
So{S, r, r') = {S + r') sign((5 + r + r') 

— TT 

dk cos{kS)ji{kr)ji{kr') = Pi{S,r,r') , 

ASr^r' 

Pi {S, r, r') = S^ {S, r,r')- S, {S, -r, r') - S, {S, r, -r') + S^ {S, -r, -r') , 



Si{S, r, r') = {S^ - 3{r^ + r'^)S - 2{r^ + r'^)) sign(á + r + r') 
With these integrations, the final expressions for Tío and Hi are 



oo />oo 



no{t) = -- dT dT' dr dr'J^{T,r)Jg{T',r')[Po{T-T',r,r')-Po{2t-T-T',r,r')] 

O J-oo J-oo JO JO 



(33) 

(34) 
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(40) 



TT r /■* /-oo ^oo , / 

^i(^) = -7^/ dT dT' dr —J4T,r)Jr+g{T',r')[P^{T-T',r,r')-P,{2t-T-T',r,r')]. (41) 

oo J-oo JO JO 



V. ASYMPTOTIC MAGNETIC HELICITY 



In the integrations for Ho(í) and T-Li{t) in Eqs. (40) and (41|, the factors Jr, Jg and J,j> have compact support. So 
the integrand is non-vanishing for some finite domain of r, r', r and r' . If we now take t to be large, while keeping r, 
r', r and r' finite, then 2í — r — r' is large and also large compared to r and r' . Then the "sign" factors in Eqs. (36) 



and (39) can be replaced by +1 and it can be checked that Po{2t — r — t' ,r,r') = and Pi(2í — r — t' ,r,r') = 0. 



Hence at late times, only the time independent pieces in "Hq and "Hi can survive, and we get 
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= dr dr' dr —J^(r,r)Jr+eÍT',r')PiÍT-T\r,r') (43) 
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nioo) = ^ V [^o(c^) + Hi(oo)] . (44) 



and 

3 V 4 

In terms of the mass of the M^-boson, mw — 3í^/2, the hehcity can be written as 

n{oo) = ^^^^^mty [Hoioo) + niioo)] , (45) 
3 

where 9^ is the weak mixing angle and for g' <C g, suiOw ~ 9^ ^ g'/g- The quantities Ho and "Hi have dimensions 
of L'^ where L is a length scale that enters the profile functions. Assuming r/ ~ r/j in the profile functions and 
left-handed hehcity (see sections below), the parametric dependence of the hehcity is 

Hioo) ^ -^^^imwrf)' . (46) 

The result can also be expressed in terms of the finc structure constant, a = e^/An « 1/137, using e = 
y/g"^ + g'^ sin 0w cos 9^, ■ 

The asymptotic hehcity depends on the time evolution chosen in ^{t) and also on the profile functions. There is no 
apparent symmetry reason for the integrals to vanish, and hence the decay of the sphaleron will lead to a magnetic 
field with non-vanishing magnetic hehcity. Below we will numerically evaluate the time-dependent magnetic hehcity 
for a choice of /i(í) and the ansatz for the sphaleron profile functions in Eqs. ([5| and (|6|. This evahiation also gives 
the asymptotic numerical vahie of the magnetic hehcity. 



VI. TIME DEPENDENT MAGNETIC HELICITY 



Eqs. pS| ) and (41) give the time dependent magnetic hehcity but the integrations are difficuh to do analytically, 
especially because the true profile functions are not known in closed form, and the time evolution function, /i(í), needs 
to be chosen. Even with the very simple profile Ansatz of Eqs. ^ and ([6|, the integrations lead to a lot of different 
terms. So we have evahiated the integrals numericaUy but this also requires a choice of the time dependent function, 
/i(í). 

One possible scheme to find a suitable /i(í) is to derive an effective equation of motion by substituting the field 
ansatz for the SU (2) sphaleron in Sec. [ll] in the energy functionaL This evaluation was done in Ref. [6 by taking /i 
to be a parameter and not a dynamical variable. The result was a "potential energy" function 

Vin) = mwiCi siv? /i + C2 sin"* /i) , (47) 

where Ci and C2 are numerical constants. If we promote /i to a dynamical variable, the energy wiU also depend on 
/i and we can obtain an cquation of motion for /i. However, this exercise is moot because it leads to a conservative 
equation of motion for /i, while the true evolution is dissipative, and once /i gets to the true vacuum value (0 or tt) it 
should stay there. One way to introduce dissipation is to replace second time derivatives by first time derivatives in 
the equation of motion. This suggests an equation of motion 

/i = Ai sin /i -f A2 sÍT? /í , (48) 

where Ai and A^ are some coefficients. The equation can be integrated and its messy solution is very similar to the 
much simpler 



Kt) = 2 



1 + tanh 



(49) 



Hence we adopt this functional form for /i(í) and proceed with our numerical evaluation of the magnetic hehcity. Note 
that /i(0) = 7r/2 and /i(oo) = tt, and hence this choice of /i(í) describes creation of baryons. 
Next we write 

nt) = ^^^^imwrf)*m , (50) 



12 3 4 

time 



FIG. 1: J-{t) versus t/rt for rt — rf = r^- 



with thc dimcnsionless function T{t) dcfincd as 

Ht)^^[no{t) + nm . (51) 

The numerically evahiated T{t) vs. t is shown in Fig. [1] for r/i = = rt. An intcresting featurc of thc plot is that 
the magnetic hehcity is negativc and conservcd at late timcs. 

Thc plot in Fig. [T] corrcsponds to Fig. 1 of Ref . [Hl where it was obtained by solving thc electroweak field equations 
on a spatial lattice starting with a perturbcd sphaleron. 



VII. SIGN OF HELICITY 



An important virtue of the present calculation is that it can rclatc the sign of the asymptotic magnetic hcUcity 
to thc dircction in which thc sphaleron dccays, which in turn is also related to whcther baryons or antibaryons are 
produced. So far we have only seen this connection by numericahy cvaluating the hcHcity integrals in the limit t — ^ oo. 
Here we wiU analyticaUy estabUsh thc conncction for smaU times, just as the sphaleron starts to dccay. 

The sphaleron is located at n — tt/2. Thcrefore we consider 



f^=l+e{t) = '^+éo5t + 0{{6t)') , 

wherc St is a smaU timc intcrval and ¿o dcnotes the velocity with which the sphaleron starts to roU. 
expand J7r, JTe and ^J^ to lowest order in St. This gives 

Jr - Fr{r)éoSt + 0{{St)^) 



(52) 



Now wc can 



(53) 



Je ^ Fg{r)éoSt + 0{iSt)' 



(54) 



Jc^ = F^{r) + OiiSt)^) , 



(55) 



wherc F^, Fg and F^ are some timc-indcpendcnt functions. Ncxt we count powers of St in the hehcity integrals. We 
get two powers from thc two time intcgrals, onc from the factors of thc currents as given above, and one power from 
the terms in square brackets in Eqs. (40) and (41 1. Thcrcforc 



nist) ^ -éoiSt)^ xl 



(56) 



8 



where T is an integral that does not depend on Since X stiU depends on the profile functions, we can only evaluate 
it numerically and, as is clear from Fig.jl] Z > 0. This tells us that an increase of \i (éo > 0), corresponding to baryon 
production, yields left-handed magnetic hehcity, and a decrease of [x (éo < 0) leads to right-handed magnetic hehcity. 
In the cosmological context, since we know that baryons prevail, this predicts that left-handed magnetic heUcity is 
prevalent in the universe [7]. 



VIII. CONCLUSIONS 



We have developed an analytical technique to study the magnetic debris produced during baryon number violating 
processes that occur via sphalerons. We have calculated the time-depcndcnt hehcity of the magnetic fields that are 
produced and found the asymptotic hehcity in terms of the decay path of the sphaleron up to quadrature. The sign 
of the hehcity depends on whether baryons or antibaryons are produced and baryon production leads to left-handed 
magnetic hehcity. Our final result for the time evolution of magnetic hehcity is shown in Fig. [1] and compares well 
with the numerical computations of Ref. [Ü]. These figures demonstrate that magnetic hehcity is conserved at late 
times. This is novel because the conservation of magnetic hehcity is usuahy discussed in the magneto-hydrodynamic 
context where a highly conducting plasma is present. In our case there is no external plasma. The (approximate) 
conservation of hehcity has also been noted in solutions to the vacuum Maxwell equations in Ref. ¡14¡ . 

The wider imphcation of our results is that there is hkely to be a unified origin of cosmic matter and cosmic 
magnetic fields. This possibility has been discussed before [31 [3 [THS] ■ The analytical technique we have developed 
makes the connection between magnetic helicity and baryon number more transparent and has the potential for 
further development. For example, one could consider an ensemble of possible decay paths described by the function 
/i(í), and the mean magnetic hehcity could then be evaluated as an average over this ensemble. Further one could 
incorporate our results in a reahstic baryogenesis scenario and derive spectral properties of the magnetic field. A 
simple analysis of this kind was done in Ref. [15], where a model for the magnetic fields produced during sphaleron 
decay was assumed. Based on our present results, it appears that the model captures most of the features of the 
magnetic fields produced in sphaleron debris, and so it is worth re-displaying the magnetic field expressions here. In 
spherical coordinates with origin at the location of the sphaleron, the model for the magnetic ficld is 



acos9 



(a2 +r2)3/2 



asiné' 2«^ - 

- + 2 („2 + ,2)5/2 (57) 

= ^e-'-/''sin0 . 
a-^ 

The extent of the magnetic field is determined by the length scale a, which is time-dependent. The simplest choice is 

a{t) = í - ío , (58) 

where ío is the time at which the sphaleron started rolhng down. As the magnetic field expands, just as in the case 
of the sphaleron, the helicity stays constant at 

Stt , , 

= - — 0.57 = -4.8 . (59) 
o 

The numerical value of the asymptotic hehcity can be adjusted by suitably rescahng the magnetic field. 

Once a space-fihing hehcal magnetic field has been produced at the baryogenesis epoch, further evolution wih 
be subject to magneto-hydrodynamical evolution in a cosmological setting e.g. [Il[l6]. We expect the presence of 
magnetic hehcity to have significant impact on the evolution. This topic, as weh as the possibihty of observation of 
early magnetic fields, deserves further investigation. 
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